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> ■ ABSTRACT 

cn ; 

^>f^ Aims. We analyze the properties of VIRGO detector with the aim of studying its ability to search for coalescing black hole binaries. We focus 
■ on the remnants of the Population III stars, which currently should be massive black holes (~ 100 - 1000 Mq), some of them bound in binary 
fN| ' systems. The coalescence of such binaries due to emission of gravitational waves may be currently observable. 
' Methods. We use a binary population synthesis to model the evolution of Population III binaries. 

Results. We calculate the signal to noise ratios of gravitational waves emitted by the system in each of the coalescence phase: inspiral, merger 
' and ringdown, and provide simple formulae for the signal to noise ratio as a function of masses of the binaries. We estimate the detection rates 
, for the VIRGO interferometer and also compare them with the estimates for the current LIGO. We show that these expected rates are similar 
^ to, or larger than the expected rates from coalescences of Population I and II compact object binaries. 
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' 1. Introduction 

, The high frequency g ravitational wave detectors are approa ching thei r specified sensitivities. The LIGO observatory 
d ' jAbramovici et alJll992 ) is already taking data, while VIRGO ( Bradaschiali 199 0) undergoes the engineering runs. It is there- 
fore very important to study in detail the detectability of different types of gravitational wave sources with the current detectors. 

A new class of sources - black hole binaries originating in population 111 stars - was recently suggested dBelczvnski et alJ 

l2004 . iBaraffe et al.' feOOl*) have shown that zero metallicity star are stable, and can form intermediate mass 100 - 500 
black holes (Heger & Wooslev 2002). They have shown that the stars with the initial masses between 140 and 260 undergo 
violent pair instability explosions and leave no remnant while the stars o utside these mass range leave black holes with essentially 
the same mass as the initial star. Numerical investigation of Bromm et al.. (.1999^ .2002 ) of Population 111 star formation indicate 
that their initial mass function allowed a much larger number of massive stars than that of Population 1. Wvithe & Loeb (20041) 
have a considered detectability of gravitational waves emitted in the process of black hole growth on the Population 111 black 
hole seeds. 

It is therefore imp ortant to ask a question whether Population 111 stars could form binaries. However, it was recently shown 
bv lSaigo et alJ ( l2004l) that a rotating metal free cloud may lead to formation of a metal free binary. Moreover, massive stars may 
form from clouds of slightly enhanced metallicity which may cool off more efficiently. All known stellar populations contain a 
large fraction of binaries. 

The detectability analysis of Population 111 black hole binaries was performed bv belczvnski et all ll2004l) for the case of the 
advanced LIGO detector. In that paper a simple model of the evolution of the Population III stars was developed. The model 
was based on numerical calculation of stellar evolution. A top heavy initial mass function was assumed, and a binary fraction of 
Population III stars of 10%. The number of stars was estimated with the assumption that a fraction of lO^^* of the baryonic mass 
in the Universe was processed in Population III binaries. 

In this paper we investigate the detectability of the Population III black hole binaries with the current configuration of LIGO, 
and with VIRGO at th e expected sensitivity. The signal to noise from a compact object coalescence for the LIGO detector has 
been calculated bv iFlanagan & Hughe sl (ll998h . In section 2 we estimate the signal to noise for VIRGO using the same formalism. 
The detection rate is calculated in section 3, while section 4 contains conclusions and discussion. 
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2. Signal to Noise Ratio 

To calculate the signal to noise ratio of gravi tational wave from binary black hole coalescence by the VIRGO detector we follow 
the formalism of lFlanagan & HughesI ill998l) . We find the average of the squared signal to noise ratio (SNR) over all orientations 
of and directions to the source, < >: 



where z is the redshift of the system, D{z) is the luminosity distance, SuiJ) is the one sided power spectral density of the 
interferometer and dEe/d fe is the locally measured wave energy spectrum. We also introduce the frequency / at the source and 
/e - (1 + zy\ f at the interferometer. In this section we follow the system of units in which G = c = 1 . 

lFlanagarr & Hughes ( 1998) approximated the sensitivity curve of the detector with the quantity /irms(./) = y/fShif), i-c root 
mean squared fluctuation in the noise at frequenc y f in a bandw idth Af - /, and then fitting them with power laws. We have 
analogously approximated the VIRGO sensitivity (|Punturol2004l) with a simple analytical model: 



oo for f < fs 

hXir' for/,</<^ 

for ^ < / < aU ^ ^ 

hU^?'^ for < f 



with parameters: 

(3) 



/s = 2 Hz 
/m = 128 Hz 
a = 3.2 
/im = 5.12-10-22 



The VIRGO design sensitivity curve and the fit ar e shown in FigureH] Our ap proximation of the VIRGO sensitivity curve differs 
from the one used for the initial LIGO detector bv lFlanagan & Hughesl(ll998h in the low frequency domain. In this domain initial 
LIGO sensitivity curve has been approximated by a function /!rms(/) ~ while the slope of the VIRGO sensitivity curve 

is steeper, with /iims(/) ~ yet it extends further to the low frequency regime. This modifies the analytical formulae for the 
signal to noise ratio, especially for the high mass merging compact objects. 

The energy spectrum d£'e[(l +z)/]/d/e has different form depending on the coalescence phase: inspiral, merger or ringdown. 
In the following we have the same definitions of the energy spectra in these phases as used bv .Flanagan & Hughes (.1998.) . For 
the inspiral phase we have used the formula: 

— = i;r2/VM2/V-'/\ (4) 
df 3 ^ ■' ^ ^ 

where fi is the reduced mass of the system, while M is its total mass, and / < /merge - 0-Q2/M, the frequency at which the inspiral 
phase ends. The merger phase energy spectrum was approximated by a flat spectrum between /merge and the quasi normal mode 
frequency /qm- = 0.13/M: 

^ = ^^^«O(/-/.erge)0(./,„r-/), (5) 
/qnr /merge 

where e„, is the fraction of mass energy radiated during merger phase which we assumed to be e„, =0.1, 0(x) is the step function, 
and FinIM) - (Ajj/M)^. Finally, the energy spectrum in the ringdown phase is 

dE A^M^f r 1 



d/ 327rV i[(/-/q„,)2 + (27rT)-2]2 
1 



[(/ + /qnr)2+(27rT)-2]2 



(6) 



^ ^-A^QM'6(f-f^,,)[\+0{\IQ)], (7) 

where Q - Trr/qm is the quality factor, t is the damping time of the quasi-normal mode of the newly formed black hole, and A is 
the initial amplitude of the wave. The quality factor is approximately Q ^ 12, we assume that the initial amplitude A - 0.4. The 
amplitude scales as F(jx/M) for non equal mass binaries. 
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Fig. 1. The VIRGO sensitivity curve expressed as ftrms(/) = ^/fShif) (dashed lines) and the approximation used in this work 
(solid lines). 



The details of the calculation are shown in the Appendix. Here we present the final results in a form useful for future 
calculations. In the inspiral phase the signal to noise ratio is: 



S M^'^ Ia[l 
— - — — \ — x< 
N Dm V M 



0.95 ^1 -3.26X 10-5M^'°/3 for < 9.9 

[.21^1-0. 



197M, 



for 9.9 <M,< 102 



5.64 X lO-'' ^8.69 x lO^^M^"'^ - 1 for 102 <M^< 2027 



(8) 



for 2027 < M, 



where = M(l + z)/Mq, and Dioo = D/lOOMpc. In the merger phase the signal to noise ratio is: 
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Fig. 2. SNR as a function of the redshifted total mass of a merging binary for VIRGO (solid lines) and initial LIGO (dashed 
lines). The dotted lines present the numerical results obtained using the approximate noise curve. The upper panel shows SNR in 
the inspiral phase, the middle panel in the merger phase, and the lower panel in the ringdown phase. 



S_ 

N 



1 4/j 
Dm M 



3.38 X IQ-^Ml'^ 



for < 9.9 



, I /m7\ Mf 4 

9.13X lO-^M.i 41og + T for9.9< 

'\ ^\9.9/ 2.0x105 3 



< 64.3 



0.25 X M, 



for 64.3 <M,< 102 



9.13 X lO-^M, ,/l +41og( - ^'^^.^/^^ for 102 < M, < 663 



(9) 



3.96 X 10"* 



for 663 <M^< 2027 



1 88 X 10^' 

9.13 X 10-^M, J- 6.25 x lO-^ for 2027 < < 13169 
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Finally in the ringdown phase the signal to noise is: 
1.80 X IQ-'^mI'^ for M- < 63.4 



S_ _ 1 An 



0.0938Mj for 63.4 < M, < 650 

4.126 X 10"^^ for 650 < M, < 13265 ^ ^ 

for 13265 < M,. 



We h ave calculated the SNR using the above approximate formulae, as well as numerically integrating the noise spectrum 
given in jPunturolbOol . The results are shown in Figure |2l We also show for comparison the SNR for the initial LIGO con- 
figuration. The inspiral phase signal to noise decreases quickly with increasing mass when (1 + z)M > lOOM© for VIRGO. 
However, for VIRGO the merger and ringdown phase S/N peaks at about 800 - IOOOMq, where the LIGO sensitivity already 
drops down. Thus VIRGO may be able to catch the signal of black hole binary coalescence with the total masses above 500 M©. 
The sensitivity of VIRGO detector is similar to that of the initial LIGO for the merging binaries with the masses below 200 Mq. 
However, the difference in sensitivity in the low frequency regime leads to improvement of the expected S/N ratio in VIRGO for 
large redshifted masses. 



3. Detection Rate 



In order to estimate the detection rate in VIRGO and LIGO we simulate the evolution of massive Population III stars, using 
the models of Belczvnski et al. (2004). We made the same assumptions about the evolution of Population III stars: the initial 
mass function is as s umed to have a slop e of -2 above lOO Mp an d extends up to 500 Mq. We use the numerical results of 
iBaraffe et all 1I2OOII) : iMarig o et al. ' fcOOl"); Heeer & Wooslev" ("2002") to describe the stellar evolution tracks, and calculate the 
He co re masses with the results of Heger & Wooslev (2002). The outcome of the nuclear evolution follows from Heger e t alJ 
(l2003h . We include pair instability supernovae so stars with initial masses between 140 and 260Mq leave no remnants. Similarly 
to the population synthesis of Population I stars, the initial orbits are drawn from a distribution flat in logarithm, and the initial 
mass ratio distribution is flat. The binary fraction is assumed to be 10%. For more details see Belczvnski et al. (2004). This 
model of evolution of Population III binaries provides an estimate of the distribution of masses and lifetimes of the binaries. We 
assume a constant star formation rate between redshift 10 and 30, and obtain the black hole binary birth distribution function 
F(M, t,z) - ®(z - 10)0(30 - z) 2 ^(M - Mi)6{t - f,), where the sum goes over all the svst ems in the simulatio n. About 10"^ of 
the baryonic mass of the Universe must have been processed through Population III stars jMadau & ReeslEoOlh which leads to 
the comoving rate of Population III star formation of R^b ~ 1.4 ■ lO^^MoMpc^^yr The calculations are carried out for a flat 
cosmological mo del with Qm - 0.3 , Qa - 0.7 and the Hubble constant //q = 67 km/s/Mpc. The luminosity distance, D{z), is 
given by (see e.g. iBulik etal]ll2004h ^ 



D(z) = c(l +z) 







df 



dz' 



dz', 



where t is the cosmic time: 



1 



//(,(! + z) VJ^m(l +z)3+Qa 
The observed black hole binaries coalescence rate is given by: 



dM " J 



d/coal(z) 1 dV 

AM 1 + z dz 



dz, 



(11) 



(12) 



(13) 



where he maximum redshift Zmax is estimated using equations (IA.2> . (IA.3> and (IA.4> for the three phases of coalescence respec- 
tively, for a given SNR value in the detector The comoving volume element is: 

dV _ AncDizf 
dz 1 + z 

The black hole binary coalescence rate, d/coai(z)/dM, depends on the birth rate of the binaries retarded by their evolutionary time: 



d/coal(z) 

AM 



F{M,t',Zf), 



where Zf is the source formation redshift, which can be found by solving the following integral: 



df 
di 



dz. 



(14) 



(15) 
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Fig. 3. Expected observed coalescence rate of Population III black hole binaries as a function of SNR for the VIRGO (solid 
lines) and initial LIGO (short dashed lines) detectors, while the long dashed lines show the expected rates in the advanced LIGO 
(.Beiczyn ski et al. 2004). The thick lines correspond to the case withour hardening of the binaries while the thin lines show the 
case when all binaries are hardened by a factor of ~ 10. Left panel corresponds to the merger phase, and the right panel to the 
ringdown phase. 

In the calculation of the binary inspiral time we also consider the possibility that there are additional factors leading to tightening 
of the orbits. These may be due to interaction with the ISM or with stars in dense centers of Galaxies, where such binaries are 
likely to sink. We model these processes by considering a model in which the orbits are additionally tightened (the binaries 
become harder) by a factor of ten. 

We present the results of the rate calculation as a function of the signal to noise in Figure |3] Each plot contains two sets of 
curves: the rate for the VIRGO detector, the rate for the initial LIGO. These rates are lower by a factor about ~ 10"* than the ones 
calculated by Belczvnski et al. (2004). In each case we present the rates assuming that the binaries evolve without interacting 
with other stars or gas, and with the assumption that they are hardened by a factor of 10. 

The detection rates do not depend strongly on the assumed evolution of binaries and are similar in the merger and the inspiral 
phase. For the assumed S /N ~ 10 the detection rate without hardening is R{> S /N) » (1 - 10) yr"'. It drops with increasing 
required SNR and for S /N x 20 it is about ten times less (/?(> S /N) ~ (0.1 - 1) yr"'). Since the required SNR in the detector 
corresponds roughly to the maximal distance out to which coalescences are detectable the rate drops with as {S /Ny^. On the 
other hand the rate increases with the inverse cube of the detector sensitivity - the minimum detectable strain amplitude. 

In Figure |3 we present the distribution of total masses of the binaries detectable in merger and ringdown phases for the 
cases with and without hardening in the VIRGO detector. The distributions are normalized to unity. These distributions show 
the most probable masses of the coalescing binaries in the framework of the given model of evolution of black hole binaries: 
without hardening and with hardening by a factor of 10. The shape of the distribution appears to depend strongly on the particular 
choice of the hardening. Therefore the optimal search strategy would be to search the entire spectrum of the masses from 100 Mq 
to 1000 Mq, as the current models can not provide more accurate predictions. The best strategy is to look for bursts of single 
frequency signals from the ringdown phase. 

4. Summary 

We have calculated the signal to noise ratio from coalescing black hole binaries in VIRGO. The low frequency sensitivity of this 
detector leads to increased sensitivity to intermediate mass black hole mergers. We find that for the merger and ringdown signals 
VIRGO has a maximum of sensitivity for systems with the redshifted masses in the range from 500 Mq and 2000 Mq. This makes 
VIRGO an especially valuable instrument in search for the signals from the intermediate mass black hole coalescences, as well 
as captures of stars by intermediate mass black holes in clusters. 

Population III stars must have included a large number of high mass stars which should have left a significant number of black 
holes with masses above 100 Mq. Some of them could have been formed binaries. Based on the evolutionary scenarios of such 
stars we calculate the properties of such binaries and estimate the rate with which they should be detectable by VIRGO. The rate 
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Fig. 4. Differential rate as a function of observed redshifted mass of black hole binary system. The solid hnes correspond to the 
case with no hardening and the dashed lines represent the case with hardening by a factor of ~ 10. The upper panel presents the 
case of detection in the merger phase, while the lower panel shows the case of the ringdown phase detection. While the binaries 
were formed at redshifts from 10 to 30, the detectable mergers happen at redshifts below unity. 



is of the order of one detection per year for the design sensitivity of initial VIRGO and similar for the current LIGO sensitivity. 
While this number is not large by itself we note that it is comparable or even larger than t he expected detection rate of the m ost 
often considered gravitational wave sources i.e. the double neutron star binaries (Belczvns ki et al .l2002HKalogera et alJ2004l) . It 
should be noted that the number of the detectable colaescences of black hole binaries originating in Population III stars increases 
quick ly with the detector sensitivity. For the Advanced LIGO the expected rates are larger by a factor of ~ 10"^ dBelczvnski et alJ 
l2004l) ." Therefore even modest improvements of the sensitivity will make the detection of these binaries more likely. 

One should be still cautious when dealing with estimates of the properties of the Population III stars. The rate we calculate 
depends on the actual shape of the initial mass function of Population III stars, and on their binar fraction. The knowledge of 
evolution of Population III stars is only known from numerical models, and hence the models evolution of population III binaries 
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is only approximate. However, these uncertainties and the possible role that Population III star could have played in the early 
Universe makes it even more interesting to look for traces of their existence, and gravitational waves searches provide such an 
opportunity. 

Therefore we consider that it is very important to include the ringdown and merger templates in the search for merging 
binaries, and to look for signals from merging binaries across the entire sensitivity range of the interferometric gravitational wave 
detectors. 
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Appendix A: General formulae for the signal to noise ratio 

In order to find the general formulae for the signal to noise ratio we insert equations J4j' © ™d j2j into equation Q and perform the 
integration. We introduce a dimesionless variable 

(1 + z) or/n, 



,/merge 

For the inspiral phase we obtain: 



(A.I) 



S_ 

N 



^r, [l2tti/3 - |2q--i/3 - iiQ.-l/3vlO/3 _ f!^(|-)ii/3] for V < 1 
7^[a'/3(2i)ii/3-a4(|.)ii/3j for a'- < v < afjf, 



for l<v<a~ (^2) 



for af„,/fs < V 

where 

2[(1 +z)M]5/3[4^/yi^] 



For the merger phase we obtain: 



S_ 

N 



V?^j|v^^] for v<l 

^y;j41nv - |(^^ for 1<v<k 

V9^ii[4 In/c] for < i' < or^ 

-^!F;„[1 +41n(if )- ^] for o-^ < v < (^.3) 

V 'r„[a'^v-''(/^ -T)] for Ko" <v<af,„lf, 

^,Jr^J[^^r~{^] for af,„lf,<v<Kaf„,lf, 

for Kaf„lfs < V 

where 
and 

6„M(I+z)2(4^/M)2 



9^, 



I0;r2Z)(z)2A2 /„,„ge('^- !)■ 

The fraction of the mass emitted in gravitational waves during the merger phase has been assumed to be e„, =0.1 iFIanagan & HugheJl998h . 
For the ringdown phase the result is 



5^ {\+zfM^A^Q{AiilMy- 
N \20;r2D(z)%,5/,[/,„/(I+z)]' 

Note that the signal to noise in the merger phase scales as e^^, and in the ringdown phase it scales linearly with the initial amplitude A. 



